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SECTION 1 
SUMMARY 


Objectives 

The contract "Study of the Liquid-Solid Transition for Materials 

Processing in Space" has two primary objectives: 

1. To analyze the behavior of dense liquids near the solidification 
point while the liquid in question is under the influence of 
magnetic fields or near-zero gravity conditions, and 

2. to do this within the framework of existing liquid state models 
and classical field theory. 

These objectives have been accomplished, and this report presents the 

results obtained thereby. 

Approach 


The approach taken to this study began by defining four tasks: 

Task 1 - Model Identification and Assessment 

Task 2 - Determination of Alternate Approach (es) as Required 

Task 3 - Incorporation of Models into Field Theory and Thermodynamic 
Property Derivation 

Task 4 - Analysis of Convection, Diffusion and R-F Field Effects. 

The models and approaches identified in Tasks 1 and 2 were combined 
formally with classical field theory in Task 3. In Task 4, the 
equations derived were used to qualitatively analyze the effects 
which external fields have on dense liquids. 

Liquid models are concerned primarily with liquid structure (thermodynamic 
properties) or liquid transport properties. The transport properties 
considered are generally diffusion coefficients, viscosity (fluidity), 
nucleation rate and crystal growth rate. Calculation of these quantities 
from first principles is the purpose of the models. However, these 
quantities by themselves do not entirely describe liquid transport, 
particularly convection. Natural convection is described by Boltzmann 
transport theory and convection in a magnetic field is described by 
magnetohydrodynamics (MHD) which is a combination of Boltzmann transport 
theory and classical field theory. Determination of the form of MHD 
theory that is required to discuss dense liquid behavior in a magnetic 
field was the function of Task 2. The models do contribute to MHD 
analysis though, by providing an equation of state for the liquid and 
by deriving parameters such as viscosity, heat capacity and magnetic 
permeability which are important constituents of the MHD equations. 
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Approach (Continued) 

Field effects on diffusion were discussed in terms of the Free Volume 
Model description of diffusion coefficient and crystal growth rate aug- 
mented by the defining equation for isothermal compressibility. This 
equation relates the change in liquid volume to the external field- 
induced pressure change in the liquid. Gravitational pressure changes 
between the earth's surface and earth orbit were derived as were magnetic 
pressure changes in the laboratory. 

Resul ts 


A thorough search of the literature identified thirty-six independent 
liquid state models plus at least six models based on the Distribution 
Function approach to liquid state theory. Most of the thirty-six models 
treat limited aspects of the liquid state. That is, they attempt to 
explain a few liquid properties (in some eases, only one) rather than 
describing the entire liquid state. Those models which do try to explain 
a large number of liquid properties are still being tested by the scientific 
community against experiment to ascertain their accuracy and the range of 
parameters to which they are applicable. 

The Free Volume Model (see Section 2.1) was selected as the best of these 
models on the basis of its simple view of liquid structure combined with 
successful prediction of liquid transport properties, and its derived 
capability for explaining thermodynamic properties. The Guggenheim Dilute 
Solution Model was selected for description of two-component ("doped") 
liquids while the Curie Law of paramagnetism was chosen to describe the 
magnetic permeability of liquids. 

"Backup" models determined to have possible application to materials 
processing problems include: 

Significant Structure Theory (1 and 2-Component Models) 

Onsager Dielectric Theory 
Kaiser & Rosenweig Ferrofluid Model 
Kuhlmann-Wilsdorf Model 
Ookawa Model 
Tong Model 

Walls and Upthegrove Model 
Shereshefsky Model 
Two Structure Model 
Widom Ntodel 
Central Atoms Model 
Kirkwood Dielectric Theory 
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Results (Continued) 

These are described in varying detail in the Task 1 report, A Summary 
of Liquid State Models for Materials Processing in Space (Reference 1 ) . 
The remaining 20 models were considered not suitable for materials 
processing problems for one or more of the following reasons: 

1. the model is grossly inaccurate when compared with experiment, 

2. the model is not yet sufficiently developed, 

3. the mathematical description of the model involves extremely large 
sums or products (-10^3 terms) or requires other calculations 
which cannot be performed analytically, 

4. the model is designed for special applications which are outside 
the scope of this study. 

Thus, several models are required to describe dense liquid phenomena, but 
these models compliment each other in creating a basis for materials 
processing in space research. It should be noted that the need for 
several models (or approximations) would have been the same had the 
Distribution Function approach (Reference 1) been chosen, since all 
approaches to describing the liquid state are imperfect at present. 

In Task 3, the five simultaneous equations of magnetohydrodynamics were 
specialized to the case of a dense liquid in the transverse field of a 
laboratory magnet. Thermodynamic properties were derived from the 
partition functions of the Free Volume Model and the Guggenheim Model, 
and the distribution of eddy currents in a conducting liquid with alter- 
nating magnetic field was calculated. Also, time-dependent equations 
for the velocity and distribution function of diffusing molecules were 
found, though solutions to these equations were not obtained. Finally 
the explicit mathematical form of the magnetic field due to a laboratory 
electromagnet was determined. This explicit representation of magnetic 
fields in terms of spatial coordinates is required to calculate vector 
derivatives of the fields (on which the magnetic forces depend). 

The magnetic field effects on liquids are really of three types; static 
and dynamic body forces and eddy current effects. Since the materials of 
interest in this contract (group III and group V elements and their 
compounds) are paramagnetic and conducting (a -lO^ ohm"^ cm-M» the 
static body forces are completely negligible compared to the dynamic 
effects. The dynamic bod^ force arises from the convective motion of 
the liquid with velocity u in the presence of a static magnetic field, 
while the eddy current force is present only^whe^ the fi^ld changes in 
time. Both these forces depend on the term 0 x B where B is the magnetic 
induction in the liquid. 

In the dynamic case, J=v x H, while in the case of eddy currents, J is 
proportional to the time derivative of B. In either case, the fields 
can in principle be arranged so that the magnetic body force opposes 
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Results (Continued) 

the convective flow caused by gravitational body forces. However, the 
actual laboratory configuration which successfully eliminates or 
drastically reduces convection may be extremely difficult to achieve 
because of the field geometry and field strength required. Furthermore, 
the fact that the region which the magnetic field occupies is limited 
in extent also limits the usefulness of this technique for materials 
processing. In earth orbiting laboratories, magnetic field strength, 
geometry and extent considerations are of course not required to 
reduce convection. 

The analysis of field effects on diffusion (within the framework of 
the Free Volume Model) yields the interesting result that application 
of a magnetic field or removal of a gravitational field will increase 
the rate of crystal growth from an undercooled melt. The conditions 
under which this occurs, and the different methods used to derive this 
result are discussed in detail in Section 2.3. Actually, by knowing 
the effects which magnetic and gravitational fields have on volume, 
viscosity and temperature, one can use several models to predict that the 
diffusion coefficient will increase with an increase in magnetic field 
or a decrease in gravitational field. The Free Volume Model then 
predicts that growth rate increases with increasing diffusion coefficient 
and decreasing Gibbs free energy. This increase in growth rate has 
been observed in the laboratory by Schieber (see Section 2.3). 

Although the equations which predict the increased growth rate are similar 
in both the magnetic and gravitational cases, this provides only a formal, 
qualitative comparison. To determine which case would predict the most 
significant increase, and to compare the magnetic case to Schieber 's 
results, quantitative calculations are required. This also involves 
improved understanding of the way free volume is distributed in liquids 
and more accurate determination of some of the Free Volume Model parameters. 

Conclusions and Recommendations 


The major conclusions to be reached from this study are: 

1. the Free Volume Model is the most practical model for use in the 
study of materials processing in space, 

2. other models (suggested in Reference 1) are required to deal with a 
wide variety of dense liquid problems, 

3. both magnetic fields and low gravity conditions will produce an 
increase in diffusion coefficient which will result in a increased 
growth rate of a crystal , 

4. laboratory magnetic fields and low-gravity conditions produce quali- 
tatively similar results (i.e., pressure terms) with respect to 
inhibition of convection and enhancement of diffusion. 
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Conclusions and Recommendations (Continued) 

5. time-varying magnetic fields induce eddy currents in the liquid which 
produce body forces and tend to disrupt convection, 

6. much quantitative and experimental comparison is needed before the 
relative merits of ground-based magnetic field techniques and space 
processing can be ascertained, 

7. the space environment provides a perfect "textbook" laboratory for 
the study of those aspects of the liquid state which are not yet 
understood, because in space the gravitational perturbation (so 
often ignored in liquid state theory) is not present. 

Thus, it is seen that the behavior of dense liquids in space and their 
behavior under the influence of magnetic fields are formally quite 
similar. But this does not indicate which of the two environments is 
most practical for actual materials processing. Determination of the 
most practical environment will require a quantitative assessment of 
field geometries and field strengths required to completely offset 
gravitational effects. This assessment should be performed for a 
number of realistic materials processing situations. 

Before calculations of this nature can be performed however, a great 
deal of liquid parameter data must be obtained for each material of 
interest. Parameters required include: electrical conductivity, 
magnetic Reynold's number, magnetic permeability, thermal diffusivity, 
viscosity, average density, characteristic convective flow velocity, 
Lennard- Jones potential parameters, hard sphere molecular diameters and 
liquid P-V-T data. Obtaining these parameters for liquid materials of 
interest will be a major task and should begin as soon as practicable. 

It is also recommended that calculations be performed utilizing the 
equations of Section 2.3 to determine quantitatively the relative effects 
of low gravity and magnetic fields on diffusion and solidification rate 
in dense liquids. These calculations should also be compared to the 
results obtained experimentally by Schieber. Finally, it is recommended 
that a number of space experiments be defined and initiated to measure 
dense liquid thermodynamic and transport properties under low-gravity 
conditions. Specifically these properties should include temperature, 
pressure, volume, internal energy, Gibbs free energy, enthalpy, entropy, 
diffusion coefficient, viscosity and solidification rate. 
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SECTION 2 

ANALYTICAL RESULTS 


2.1 A REVIEW OF LIQUID STATE MODELS 

There are two analytical approaches aval Table for the study of the liquid 
state: the Distribution Function approach and the Liquid Model approach. 

The Distribution Function approach is the more rigorous of the two, but 
often its mathematical formulation involves integral equations or other 
computations which cannot be solved in closed form. Even when solution 
is possible, the results do not compare more favorably with experiment 
than do results obtained from the model approach. The Distribution 
Function approach is particularly inaccurate in its description of 
liquids near the liquid-solid transition. On the other hand, liquid 
models are available which are specifically designed to describe dense 
liquids near the liquid-solid transition, and which are constructed to 
be mathematically more practical than the Distribution Function approach. 
In addition, most models are based on combinations of experimental 
observations and simple assumptions which make them conceptually easy to 
understand. For these reasons, the Liquid Model approach is recommended 
for application to practical materials processing problems in the space 
environment. 


In general, different models describe different liquid properties. For 
instance, some models describe thermodynamic properties while others 
describe transport or electric or magnetic properties. A few models 
such as the Free Volume Model (as developed by Boeing) or the Significant 
Structures Model, address both thermodynamic and transport properties. A 
detailed discussion of liquid state theory and liquid models is contained 
in the Interim Report for Task 1 of this contract (Reference 1). However, 
a brief review of liquid state concepts, assumptions and models important 
to the study described by this report is appropriate here. 

Basic Assumptions 

The most important assumption in the study of liquids is that almost all 
the properties of a liquid may be calculated from the interactions of 
pairs of molecules, and that contributions from "3-body" or higher order 
interactions can be neglected. Equally important is the assumption that 
these interactions are additive. Specifically, it is assumed that the 
total potential energy of a system of molecules is given by: 


$ 


N N 


*j n 

^ 1=1 j=l 


( 1 ) 
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2.1 (Continued) 

where tt>ij (rij)is the potential energy between molecules 1 and j, which 
are separated by a distance While it is obvious that the force on 

a molecule in a liquid is due to the effects of all other molecules in 
the liquid, the pair approximation (equation 1) is necessary to simplify 
mathematical analysis and has been justified by agreement with experiment 
in many instances. 

Thermodynamic and mechanical properties of liquids are derived from a 
quantity called the partition function (P.F.) which has the form: 

2 = Z|^Q(N, V, T) . (2) 

The first factor, z., is the kinetic P.F. which is derived from the 
motions of the individual particles in the liquid. If only translational 
motion of the particles is considered, the kinetic P.F. is simply: 

, 3N 

^ _ 2TrmkT I 2 _ -3N ... 


where T is absolute temperature, N is the number of particles, k is 
Boltzmann's constant, h is Planck's constant, m is the mass of one 
particle and a is the thermal wavelength. 

The second factor in equation 2 is called the configuration partition 
function and is defined by the expression: 


Q(N,V,T) = Jr, / e d3x (4) 

n. y 

which is an integration over the three space coordinates of the N particles 
in the system where 4' is the total energy of the system. As equation 1 
indicates, for real liquids 4 is an extremely complicated function of 
coordinates. This is why many different models of the liquid state have 
been developed - to allow one to calculate approximate expressions for 
Q(N,V,T) from certain simplifying assumptions about the liquid situation 
of interest, which will then yield predictions about that situation which 
are reasonably accurate. 

An assumption usually made to evaluate equation 4 is that all tf). . have 
the same form. Then equation 1 reduces to: ^ 


N 

4= E <^(r.) (5) 

i=l ^ 
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2.1 (Continued) 

where i now labels the distance vector, r., between each pair of molecules. 
In cell models of liquids it is assumed that all molecules are situated 
in cells defined by their nearest-neighbors. The average radii of all 
cells are assumed to be equal, that is r. = r., so that (Reference 2): 

* J 


$ = N(j)(r) 


( 6 ) 


This greatly simplifies the integration over coordinates and thus the 
configuration P.F. itself . Unfortunately, the correct form for 4> (r) 
is not known, though several approximate forms which are based on careful 
experimental observation are available, The most widely used of these is 
the Lennard-Jones potential: 


r . 12 

6 7 

(2.1 

/ a \ 

1 r 1 

1 r ] 


where eis an interaction energy constant and is a distance parameter 
equal to the value of the intermolecular separation, r, which makes 
<j> (r) = 0. Thus the assumptions of equations 5, 6 and 7, coupled with 
the specific assumptions of a given model, allow one to obtain approxi- 
mations for Q{N,V,T). 

Finally, the thermodynamic properties of the liquid are obtained from 
the partition function by the following set of equations: 


FREE ENERGY: 

F - - 

kT In Z 

( 8 ) 

INTERNAL ENERGY: 

— - 1 
kT ' 

- 9lnZ' 

r^T V 

( 9 ) 

EQUATION OF STATE: 

It 

lalnZ] 

1 1 

T 

( 10 ) 

INTERNAL ENTROPY: 

|=mz+T 

|9lnZl 

1 ) V 

01 ) 

ISOTHERMAL COMPRESSIBILITY: 

3 = 

. 1 (kT r’ 

V [ ' 9V2 Ij 

02 ) 
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2.1 (Continued) 

Once these basic functions are known, all other thermodynamic properties 
may easily be found (References 3 and 4). The next three subsections 
will discuss the three models most applicable to the Liquid-Solid Transi- 
tion Study. 

The Free Volume Model 


The Free Volume Model, which has been developed by Turnbull and Cohen 
(References 2 and 5-7), is perhaps the most simple model of the liquid 
state which is at the same time realistic in concept, practical in its 
mathematical formalism and reasonably accurate in its predictions of 
liquid properties. Development of this model has been motivated by 
observed liquid properties whi chi relative to solid properties, appear 
to be strongly correlated with the corresponding volume differences. The 
Free Volume Model makes no specific assumption about molecular shapes, and 
thus leaves one free to assume spherical shapes, or more realistic shapes 
when considering polyatomic molecules, depending on the accuracy one 
wishes to achieve. 

The basic assumption of this model is that a liquid molecule is caged 
most of the time by its nearest neighbors, and that within this cage 
its free volume is defined by: 


Vf = V -Vq (13) 

where v is the specific volume (or volume of the cage), V/N, and v is 
the actual volume of the molecule, v is dependent on pressure, though 
for normal pressures this dependence °is negligible. It is also assumed 
that the potential energy of the molecule in its cage is approximately 
equal to the intermolecular potential, <i)(r), which is chosen to be the 
Lennard- Jones potential, equation 7. The total partition function for 
this model is then: 


Z(N,V,T) = 


_- 4 .(r)/kT 
X3 ® 


(14) 


If the procedure of equation 10 is used to derive the equation of state, 
one finds for the Free Volume Model: 


PV = 


NkT MW M 

I-Vq/V 9V 


( 15 ) 
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2.1 (Continued) 

In order to evaluate the partial derivative, we first note that: 


NV 


9(j) 

av|T 



( 16 ) 


Further work with the Free Volume Model at Boeing has centered around the 
"lumped-parameter" assumption that the specific volume is directly pro- 
portional to the cube of the intermolecular separation with proportion- 
ality factor a : 


V = a^r^ . , (17) 

This allows complete evaluation of the partial derivative in equation 15 
to yield: 


PV 

- V ^ +1 

6 

1 3a] 2 [bo U 


|3al 

j 2 j 


[2 1 

Nk'T 

~ 1-Vo/V 

^ kT 

[ Ztt 1 \ V / 


j 


iVi 

1 ^ 2 


where a = a /N (N is Avogadro's number) and b = (2/3)'it N If 
equation 18°is°solved for a^, then: ° 


a 


2 = 




•si 

e 


I-Vq/v 


-W 


(19) 


where W is the compressibility factor, PV/NkT. Thus a depends on temper- 
ature and pressure, but because of the form of this dependence, the varia- 
tion of a from a constant value (for a given liquid) over certain ranges 
of T and P should be small. Table I gives a values and deviations for 
ten liquids over temperature and pressure ranges corresponding to the 
dense liquid phase. 

Equation 18, then, is the equation of state for the Free Volume Model, 
and its accuracy is determined largely by the deviation of the a 
parameter over the range of T and P to which it is applied. Table I 
shows that this deviation is less than 20% for all liquids studied and 
is less than 1% for water. A 20% inaccuracy, while not desirable, is 
still normal for theoretical predictions of dense liquid properties. 
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2.1 (Continued) 

Since the derivative of intermolecular potential with respect to volume 
occurs often in calculating thermodynamic properties from the partition 
function, relation (17) has allowed the derivation of 


INTERNAL ENERGY: 


E 3 . Ae 3a\ 
nFT" 2" FT 27^ ,v , 


3a 


2ir 


"• b, 


' 1-1 


( 20 ) 


CHEMICAL POTENTIAL: 


kT 


^ = In ^ + 


1 + 4e f 3a 


Vf 1 -Vq/v kT |27t 


2 I . V2 
V \ 


3a 


L I 


bol2 


-3 


1 ^ 4e/ 3a] 

2 

bo' 

2 

r j 3a °o ~ 

1-Vq/v kT\2TrJ 


ivJ 


( 27/ ( V / '^J 


FREE ENERGY: 
ENTROPY: 


(3al"[bol'f 

i 3aT 

) bo 1 

2 

-1 

\2-.j Iv 1 _ 

2ir 

/ 

1 V i 


( 21 ) 
( 22 ) 

(23) 

(24) 


The Free Volume Model treats diffusion by assuming that random density 
fluctuations can open holes in the cages of individual molecules which 
are large enough to permit the translation of a molecule across this void 
created within its cage. Thus diffusion is viewed as the result of a 
redistribution of free volume within the liquid rather than as the 
result of an activation in the ordinary sense (Reference 5). If the 
free volume exceeds some critical value, v*, then diffusion can take 
place, the diffusion coefficient being: 



V* + ^ e"’^'^*^''f 

Y i 


(25) 


where 6 is the ratio of the mean-free-path to the "instantaneous" free 
volume, u is the velocity of the molecule as it diffuses (taken to be 
the gas kinetic velocity) and ^ is an overlap factor, lying between 1/2 
and 1, which accounts for the fact that more than one molecule can have 
the same value of the free volume. 


Once the diffusion coefficient is known, several other transport pro- 
perties become available. First, the viscosity is inversely related to 
the diffusion coefficient by: 
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2.1 (Continued) 

- kT 

3ttR^D 


( 26 ) 


where R is the molecular diameter. Crystallization is characterized 
by two fate functions: the frequency per unit volume at which crystals 

nucleate, I, and the crystal growth velocity, u . These functions are: 


I 


nd'^ 


e' 3F2kT 


u 


c 



e 


AG 

kT 1 


(27) 


(28) 


where: oi is a parameter identified with the liquid-crystal surface 

tension, 

F is the free energy of crystallization per volume, 

a 6 is the change in free energy at the crystal -melt surface, 

and f is the fraction of lattice sites in the boundary to which 
° molecules can be attached. 


Strictly speaking, D' and D" are rate constants from nucleation theory, 
but Turnbull [Reference 2) finds that to a good approximation, D", D' and 
D can be considered identical. 

It should be noted that definitions of the free volume other than equation 
13 do exist (References 2, 7 and 11), However, they all involve some 
combination of the specific volume, which is a thermodynamic variable, 
and some reference volume which is approximately constant. Therefore, 
since v can be calculated from handbook data for various molecules, 
equation 13 is probably more practical than the other definitions. 

The critical volume for diffusion, v*, is not really specified by the 
theory. But it is found experimentally that the product, yv*, is near 
the molecular volume, v , so that calculations of transport and kinetic 
properties in the Free volume Model are greatly simplified. 

T he Guggenheim Dilute Solution Model 

In his classic book (Reference 12), Guggenheim presents a general treat- 
ment of the dilute solution (or "doped" liquid) situation which leads to 
a simple partition function for the mixture when the P.F. for the solvent 
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is known. The initial assumption is that interactions between solute 
molecules can be neglected, so that these molecules are viewed as a 
perfect quasi-gas moving freely in a region of constant potential 
energy, x . x depends on the nature of the solute, the nature of the 
solvent, ^and ^the temperature. The partition function is thus: 


Z 






(29) 


where N is the number of solute molecules, z,^ is the kinetic P.F. for 
the solute (including internal effects, as desTred), and Z° is the P.F. 
for the one-component solvent. V is the total volume, as usual. From 
equation 29, the free energy is simply: 


F 


-k.T In - NgXg 


NgkT In 


Ns 


(30) 


Differentiating F with respect to N yields the chemical potential for 
the solute: ■ 


y 


s 


kT 



In 




(31) 


where N, is the number of solvent molecules and V, is the volume of the 
solvent. By using the Gibbs -Duhem relation (Reference 13), the chemical 
potential for the solvent is: 


Pi. 

kT 


kT 


+ In 


U. 


N, + N, 


(32) 


where the subscript 1 refers to the solvent and the superscript o indicates 
a quantity for the pure solvent substance. The importance of the chemical 
potentials arises from the fact that in mixtures, the chemical potentials 
are identical to the partial molar Gibbs potentials of the components 


Pi = Gi , 


( 33 ) 


This is useful for two reasons: 

(1) in any reversible isothermal process (such as the liquid - solid 
transition) the decrease in Gibbs potential is equal to the net 
work done by the system; and 
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(2) when a mixture undergoes a first order phase transition, the 

equilibrium distribution of the substance s between the two phases 
is determined by the condition that pg have the same value in each 
phase, since equal Gibbs potentials is the criterion for phase 
equilibrium. 

Guggenheim considered only thermodynamic properties of dilute mixtures in 
his model, and neglected the transport properties. However, several other 
models (Reference 1 ) , including the Free Volume Model, provide descriptions 
of transport properties of multi -component mixtures. Thus other models 
must be used in conjunction with the Guggenheim model to completely des- 
cribe the liquid system. The thermodynamic functions for dilute solutions 
which have been derived from the Guggenheim model (expressed in terms of 
the pure solvent functions) are: ■ 


FREE ENERGY: 

EQ. OF STATE: 
ENTROPY: 

INTERNAL ENERGY: 
ENTHALPY : 
SPECIFIC HEAT: 


F 

_ F° 

NikT 

NikT 

PV 

pov 

NikT 

NikT 

S _ 

: J1 + 

Nik “ 

Nik 

E 

E° 

NikT 

N^kT 

H 

HO 

NikT 

NikT 


Ng rx 


^ + In 




3 + In Ve 

2 "'"^3 


2 kT 


[ 2 kT J 


C = c° + - N k 
V V 2 s 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


The Curie Law Of Paramagnetism 

The Curie law of paramagnetism was discovered experimentally by Pierre 
Curie in 1895 and derived theoretically by Paul Langvin in the early 
twentieth century. Langevin's procedure was to derive the average value 
of the molecular magnetic dipole moments in the direction of an applied 
field. The result is: 

m^H. 

m cos e = (40) 
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where H. Is the resulting internal field, 0 is the angle between the 
field ana an individual dipole, and m is the permanent dipole moment of 
a molecule. 

The volume magnetization of the bulk material is obtained by multiplying 
equation 40 by the density, n. The magnetic susceptability of the material 
is then: 


- JL - 111 

‘ " 3kT 


(41) 


Once the susceptibility of a liquid is known, the magnetic permeability is 
given by: 


w " X|^) 1^0 


(42) 


where uo is the permeability of free space. The magnetic induction, B, 
which determines the force on a material in an external field, “fi, is 
then related to the field by: 


B = p H 


(43) 


The Curie Law model is important in materials processing because all 
natural liquids, that is those which are not suspensions of ferromagnetic 
powders, are paramagnetic or diamagnetic, with most materials of interest 
in space processing being paramagnetic. This means that equation 43 holds 
whereas it is invalid for ferromagnetic material. Equation 43 in turn 
allows straightforward calculation of the effect^external magnetic fields 
have on liqu4.ds. This would not be possible if B were not a linear 
function of H. In practical calculations, it is not always necessary to 
make use of equation 41. Since is on the order of 10"® for para- 
magnetic materials (-10"® for diamagnetic materials), it is often 
sufficient to use a constant value for u, since it will change only 
slightly with varying temperature and density. In any event, it is the 
magnetic permeability of a liquid which couples the liquid directly to 
a static magnetic field and, with the electrical conductivity, helps 
couple liquid and field when motion is involved. 
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2.2 FIELD EFFECTS ON CONVECTION 

To understand the effect external fields have on convection, one must 
first_ understand the general theory of fluid flow, which is derived from 
Kinetic Theory (Reference 14). The first step is the derivation of the 
Boltzmann transport equation which is a n§nl;i^near integro-differential 
equation for the distribution function f(r, u, t)* f is defined so that 


f (r,u,t) d^r d^u 


is the number of molecules at time t whic^ have positions lying within 
a volume element djr in the region about r and within a velocity-space 
element d^u about u. Needless to say, solution of the Boltzmann trans- 
port equation poses a tremendous mathematical problem. However, by 
combining approximations to the transport equation with the conserva- 
tion theorems for mass, momentum and energy, the equations of viscous 
fluid dynamics may be derived. These are commonly called the continuity, 
Navier-Stokes and temperature equations and are usually combined with 
an equation of state which provides an independent relationship between 
temperature, pressure and density. 

It should be noted here that density, rather than specific volume, is 
customarily used in discussing fluid flow problems. These variables 
are related by: 


n 




(44) 


where N is the total number of molecules and V is the total volume of 
the system, n is sometimes called the number density to distinguish it 
from the mass density, p, and m is the molecular mass. 

Boltzmann transport theory and the fluid dynamics equations derived from 
it might be considered as an additional model to those in Section 2.1. 

Or they might be considered an alternate approach to the problem of 
field effects on liquids in the sense that fluid dynamics is quite 
different from the models in Section 2.1. In this sense, transport 
theory may be considered as part of Task 2 of the contract, which was 
to "determine alternate approach (es) as required." 

The Equations of Maqnetohydrodynamics 

If a liquid is placed in a magnetic field, the theory which describes the 
liquid-field system is called magnetohydrodynamics (MHD). For the descrip- 
tion of such a system, one adds to the four fluid dynamics equations 
discussed above (modified to include magnetic forces) a fifth equation 
which governs the behavior of the field in the region of the liquid. 

The five simultaneous MHD equations are then: 
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CONTINUITY: |^ + v-nu ^ 0 (45) 

P - ■^•uj+ n72u (46) 
TEMPERATURE: + u-^ T ^ 4 (v-u) T + aV^T + n$ (47) 

Ou I Cy n 

STATE: P = T^-fT " + Agn^ (48) 

I "VQn 

FIELD: L |H = ,2j5 4 rJ, j M j (49) 


where: u is the velocity of a small liquid volume element, 

fj is the total body force (force per unit volume) due to 
external fields, 

P is pressure 

n is viscosity in poise (1 poise = 10^ ) 

T is absolute temperature 

Cy is specific heat at constant volume 

IC 

a is thermal diffusivity a - -;r- 

pCy 

n4>^ is a viscous heat dissipation term 
Vq is the hard-sphere volume of a molecule Vq - 
V is magnetic diffusivity ^ 

Rjj^ is magnetic Reynolds number (Rj^ = Luya) 


A 3 and A 5 are constants related to the Free Volume model parameters 
(Section 2.1 ) by: 
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Thus, the models discussed in Section 2.1 enter the MHD equations at 
three points. Of course equation 48 ccrnies from the Free Volume Model 
(it should be replaced by the equation of state from the Guggenheim 
model, when dealing with two-component liquids). The magnetic 
permeability occurs in both v and Rm in the field equation (49), and 
also contributes to the external force term, f-|-, in the Navier-Stokes 
equation. This body force term will now be discussed in some detail . 

Explicit Terms and Approximations 

In the case of interest, i.e., a liquid under the influence of both 
gravitational and magnetic fields, the force term can be considered 
the sum of three component terms: 


GRAVITATIONAL: 
STATIC MAGNETIC: 
MAGNETODYNAMIC: 



(51) 

(52) 

(53) 


f, is obviously proportional to the local gravitational acceleration, 
which can be expressed as the gradient of a potential function, ^ (or 
of pip where density gradients are small). The static term, fg, is the 
body force which is applied to the liquid when the field is constant in 
time and when the liquid is at rest. Since y has the form of equation 42, 
the derivatives of y can be considered to be derivatives of the suscep- 
tibility, xm* paramagnetic (or diamagnetic) liquids, xm is on the 
order of ±lu"^ and the derivatives of xm correspondingly small. So 
to a good approximation we have: 



(54) 


The magnetodynamic term really applies to two different cases depending 
on whether the magnetic field is constant in time or not. For the case 
of a static field with convection in the liquid. Maxwell's fourth 
equation states: 


J = V X H 


(55) 
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Therefore equation 53 becomes: 


fs = p(vxH)xH = y[(S-v) H - ^ V H2] , (56) 


by applying a vector identity (Reference 15). We can now write the 
entire external body force term as: 

fj = y(H-v) H - V [pi(i + ^ y H^] . (57) 

In simple geometries, such as those in which H| is not a function of x^ , 
the first term of 57 vanishes. This allows a clear comparison of the 
gravitational and magnetic forces. The resulting expression; 


fy =5 - V[pT() + ^ y H2] (58) 

says that p^^;, which can be called a gravity-induced pressure, and 
which is a magnetic pressure, are completely equivalent in a format sense. 
That is, magnetic fields produce body forces the same as gravity. This 
means that if the magnetic field can be controlled so that: 

-V (^ yH^) = + V (pill) (59) 


then the external body force can be minimized. The pressure gradient 
in equation 46 can also be interpreted as a body force, and the magnetic 
field can be used to balance this force as well. Conversely, if the 
liquid is in a low-gravity orbiting laboratory, and a net body force is 
desired, a magnetic field may be used to replace gravity by equation 58. 

Of course, for a rigorous numerical solution of MHD problems, one must 
solve equations 45-49 simultaneously for the density, temperature and 
velocity distributions and for the equilibrium pressure and magnetie 
field values subject to the boundary conditions of the specific problem. 
However, the equations of state and continuity are basically constraint 
equations (if pressure is approximately constant). Therefore, qualitative 
information about fluid flow patterns, temperature distributions and 
magnetic field behavior can be obtained from the remaining equations if 
it is assumed that all other quantities in the respective equations 
(46, 47, 49) are known. This will be attempted in the following para- 
graphs. 
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Magnetic Field Representations 

Before considering how magnetic fields affect liquids in more detail, it 
will be helpful to digress into how laboratory magnetic fields are 
represented mathematically. 



FIGURE 1 : MAGNET GEOMETRY 


Figure 1 depicts the two pole faces of a laboratory electromagnet. The 
equation for the magnetic induction between the faces (before a sample 
is introduced) is: 


v2 =0 . (60) 

The general solution in cylindrical geometry is: 

^ lo (ikr) (61) 

where k and Ck are constants determined by the boundary conditions, r is 
the radial distance from the central axis, and Iq is the zeroth order 
Bessel function. 

If the magnetic surface charge, M, is now calculated from the currents 
in the pole faces: 


M 



V • (r X Jo) 


(62) 
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2.2 (Continued) 

equation 61 may be expressed as an expansion in either r or z (Reference 
16): 

= 4ttM[Fq + F 2 (z/d )2 + F 4 ( 2 /d)‘' +...], z2 < d2 (63a) 
= 4 ttM[Fq - l/ 2 F 2 (r/d )2 + 3 / 8 F 4 (r/d)^ - . . .],r < d (63b) 
where the coefficients 

Fq = 1 - COS0O 

¥2 = 3/2(cos36q - cos^0q) (63c) 

F 4 = -5/8(3cos^0q - lOcos'^ 0 Q + 7 cos^ 0 q) 


The explicit dependence of Bz on r and z shown here allows the calculation 
of vector derivatives of the field (on which the forces depend) if one 
recalls that Bz = liHz. 


Magnetic Diffusion 

Equation 49, which is derived in Appendix A, relates the spatial dependence 
of the magnetic field in a moving fluid to the time dependence of the field. 
Problems involving dense liquids can be greatly simplified by considering 
the magnetic Reynolds number in this case, where, from equation A- 8 : 


Rffi — Lu |io . 


For typical group III- group V elements, the properties y and a have the 
following orders of magnitude: 


y “ 10 "®f^sec/cm 


0 = lO^ f2"^cm"^ 

If the experimental liquid sample has dimensions on the order of 10 cm 
and natural convection in the liquid produces a velocity on the order of 
10 "^ cm/sec, we have: 
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Rpi = (10cm)(10"^cm/sec)(10“®J^sec/cm)(10%"^cm"^) = 10“^. (64) 

Since Rj,, is so small, the second term on the right hand side of equation 
49 can be neglected. This equation then becomes: 


i ^ 

V at 


V2 


(65) 


which is a simple diffusion equation for the way the field penetrates 
the liquid. The time constant for this penetration will be: 


T = ycL^ 



( 66 ) 


Using the previously assumed values, this is: 


T “ lO'^sec = 1 msec. 


Thus for static fields or for oscillating fields with frequencies of 
less than 10^ Hz there will be negligible screening of the field. How- 
ever, for higher frequency fields, for instance RF fields at 10® Hz, the 
fields will be screened with a penetration or "skin depth" of: 



For the above values of y and o, RF fields have a skin depth on the order 
of 0.4 cm. This is the region in which eddy currents will be of greatest 
magnitude for such fields. 


Eddy Currents 

The screening, or exclusion, of external magnetic fields in a conducting 
medium is accomplished by the electric currents flowing in the region 
of the medium (liquid) defined by the skin depth. For a conducting liquid 
in a cylindrical container placed in a transverse magnetic field (see 
Figure 2), there will be no ordinary current flow, since the liquid is 
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not part of an electric circuit. However, if the magnetic field changes 
in time with a frequency w, eddy currents will be induced which in turn 
will give rise to a magnetodynamic body force: 

^3 = Je* X t . (53) 

Here ?e* ''s the complex conjugate of the induced eddy current, Je* 


z 

I 




FIGURE 2: LIQUID CONTAINER-MAGNET GEOMETRY 


Eddy currents are determined in the following manner. According to 
Faraday's law, a changing magnetic field gives rise to an induced electro- 
motive force, Eq, which is related to the field and the eddy current by: 

Ee = 1 § Je • dt = -//s || • ds . (68) 

This states that the line integral of the current around its circuit 
loop is equal to the surface integral of the rate of magnetic flux change 
through the surface enclosed by the loop. Equation 68 also assumes that 
Ohm's law: 



( 69 ) 
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Is valid for the liquid material under consideration. The eddy current 
must also of course satisfy Maxwell's fourth equation which is repeated 
here: 


= ^ X 


(55) 


From these four equations (53, 68, 69 and 55), the shape and direction 
of the eddy current loops and the liquid flow patterns they cause will 
now be deduced. 

First of all, for the geometry assumed in Figure 2, the magnetic field 
will have the form: 


= [Hr(r,e,z) f + H0(r,e,z) e] e^“^ 
which yields, from equation 55, the current expression: 

3e ^ 


- 0 4 - 

(aHe 

1 \ . 

dz 3Z 

\ 8r 

■ r 1 ^ 


.iut 


(70) 


(71) 


Equation 71 indicates that in the geometry of the liquid, the currents 
will have no simple symmetry. But from equation 68, we know that the 
current loops will lie in planes normal to the field direction (see 
Figure 3). Note that the direction of current flow will reverse when 
the oscillating field reverses its direction. Of course, the currents 
are much stronger in the skin depth region than they are outside this 
region, i.e., for r < (rm-6) there is practically no current. Thus, with 
the exception of those small loops which fit entirely in the skin depth 
region and are nearest the pole faces, the current is strongest in the 
±z directions. 
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The ±z components of the eddy current loops, then, produce the largest 
body forces. In the absence of all other external forces, equation 53 
tells us that these eddy current-produced body forces will cause liquid 
flow patterns to be set up in the r-e planes over the length of the 
liquid cylinder. It should be noted that the body forces produced by 
eddy currents are independent of time since the imaginary exponential 
time dependence vanishes upon performing the operation x t. This 
means physically that the field and current oscillate in phase to give 
a constant force. 

Before determining liquid flow patterns, we first need to know what the 
skin depth configuration is. From the surface integral in equation 68, 
it is^obvious that the net magnetic flux inducing the eddy currents is 
§ • ds - B • rds = Bds cose, since r is the unit vector normal to 
the cylindrical surface. The greatest contribution to the integral will 
then be near 0=0, that is, in the direction parallel to the field, while 
the contribution will fall off toward the direction transverse to the 
field, 0=tt/2. Thus the skin depth region will assume the shape shown 
in Figure 4 where the liquid cylinder is viewed from the top. 


8 = ir/2 



FIGURE 4: SKIN DEPTH CONFIGURATION AND- LIQUID 

FLOW PATTERNS 


From the current and field directions shown in Figure 3, equation 53 indi- 
cates the force (and therefore the liquid flow if no other forces are pre- 
sent) direction to be perpendicular to and lying in the r-0 plane. Since 
the force will be most intense in the skin depth region, the net effect 
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is to produce a flow going from the 0 = ± tt /2 poles toward the 0=0 or 0=11 
points, thus yielding the patterns illustrated in Figure 4. These 
patterns are repeated in each r-0 plane of the liquid cylinder. 


In physical situations, there are always forces other than magnetodynamic 
forces present. Specifically, in a gravitational field, buoyancy-driven 
natural convection currents will be flowing. Whatever the shape of the 
convection flow patterns, they will have a large ±z component which will 
be disrupted by the eddy current-produced flow described above. In 
addition to breaking up convective flow, eddy currents produce heating 
in the region they inhabit with the heating rate having the proportionality 


Q( watts) oc F] 



(72) 


where Fi(rm/6) is a function which depends on the geometry of the system 
(Reference 18). This heating produces a temperature rise in the liquid 
which has an effect to be discussed in Section 2.3. 


Magnetic Viscosity 

By using equations A-3 and A-4 from Appendix A, it can be shown that: 

J X t = ya (tt X - yH^u^) (73) 

where i^t s'"® the electric and velocity field components transverse 

to the magnetic field. The Navier-Stokes equation l46l can be written as : 

+ P0(tt X B - (74) 

where du/dt is the total time derivative of velocity, and 


z = -V (pip + P -[ri/3jv • u) + nV^u . (75) 

In liquids it is possible for z to be small, and when eddy currents are 
the largest currents present, tt will be negligible. Under these con- 
ditions, induction drag will cause the motion across the field lines to 
decay with a time constant: 


X 



(76) 


For a system of density 5 gm/cm^ with a = 10^ Q“^cm“^ and B = 10^ gauss, 
the time constant is on the order of five seconds. This drag can be 
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interpreted as being due to a magnetic viscosity with an effective force 
comparable to: 

Fn, = ay2H2u . (77) 

This corresponds to an effective magnetic viscosity: 

= ctu2h2l2 = (78) 

which, for our example a, L and B is roughly one poise. For comparison, 
the viscosity of mercury is 10“^ poise while molten glass at Its working 
temperature is 10^ poise (Reference 2), The Hartmann number is defined 
as : 


h = Vt 

and is the parameter used to determine the Importance of magnetic vis- 
cosity which dominates for h > 1. 


For the case i§ which tt '■s not negligible, the lines of force (always 
transverse to H) will be frozen into the liquid and parried along with 
it. If the velocity of these force lines is called w, we have for the 
current in the moving system: 




so that: w 


0 (?t + y w X ff) 
Et X 


(80) 


This is the so-called "E cross B drift" which, when substituted into 
equation 74, shows that ^ 0, the liquid motion approaches in a time 
T a velocity such that ut = w. Hence the transverse motion is still 
destroyed unless ut/L >> 1. 

If ? is not negligible, then relative transverse motion is not destroyed 
but approaches a value ut such that: 

\ = ay2H2 {tit - W) . (81) 

The force ?+ tends to push the lines of force past the field lines at such 
a velocity that the induced drag (viscosity) just balances i. Thus there 
is always adjustment of the liquid motion (in a very short time) In such a 
way as to produce a balance of the viscous forces I The result is a 
steady-state Newtonian- type flow (after time t > t) with: 
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[4 (p^ + p -{n/3]v • u)]^ + n(v2j)^ = Gp 2 H 2 (^^ _ _ (82) 


Recapitulation 

The point of these discussions of magnetic diffusion, eddy currents and 
magnetic viscosity is to provide a theoretical explanation of the effects 
which magnetic fields hav§ on dense liquids and how they relate to the 
internal driving force, -vP, and the gravitational body force, pg = -Vp 4‘. 
It is to be remembered that the major effects in counteracting gravity- 
induced convection are: 


(1) eddy-current produced body forces which tend to disrupt convection 
flows, and 

(2) a magnetic viscosity (equation 78) which increases viscous drag 

throughout the liquid and rapidly decays all motion transverse to 
the field lines for situations in which x and small. 
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2.3 FIELD EFFECTS ON DIFFUSION AND SOLIDIFICATION RATE 

Diffusion Coefficient Increase 


There are many liquid state models which derive expressions for diffusion 
coefficients from first principles (References 2 and 19). The major 
diffusion coefficient models fall into three categories, and these are the 
ones which are relevant to the present study. The first category contains 
the Free Volume Model in which the diffusion eoefficient is primarily a 
function of free volume: 


D = M (yv* + v^) (25) 

where as before S in this expression is the ratio of molecular mean-free- 
path to the "instantaneous" free volume. In the second category are 
models which relate diffusion coefficient to the ratio of temperature to 
viscosity, as in the Stokes-Einstein model: 

D = k2 ^ . (83) 

In the third category, represented by Swalin's model (Reference 20), the 
diffusion coefficient is a function of T^ only: 


D = k3 T^ . 


(84) 


Based on the findings of Section 2.2, it is easy to derive the effect 
magnetic fields have on the diffusion coefficients in these last two 
categories. In Section 2.2 it was shown that both the temperature and 
viscosity of a conducting or semi -conducting liquid will increase when 
a magnetic field is applied. The result for the Swalin model is 
trivial. If D' is the diffusion coefficient in changed field conditions, 
it is obvious from equation 84 that D' > D. 

Things are not as simple with the Stokes-Einstein type models, however. 

In this case, the field-induced diffusion coefficient may be written: 


D' = k 


T + AT 


2 r) + An 


(85) 


where aT and An are the field-induced terms. Therefore we have: 


D' _ 1 + AT/T 
D " 1“ + An/n 


( 86 ) 
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Hence if aT/T > An/n, D* will be greater than D. However, it is not obvious 
that this is always the case. Of course, An/n is simply the square of the 
Hartmann number, but to date, no straightforward way of calculating aT/T 
has been developed. 

Now consider the Free Volume Model. Here we have: 

y 

D' = (yV* + v^ + AV^) e"Vf+AVf (87) 


and a little observation shows that D' is greater than D if one assumes 
that application of a magnetic field causes an increase in free volume. 
This will presently be shown to indeed be the case. At any rate, both 
the Free Volume Model and Swalin's model indicate that magnetic fields 
definitely increase the diffusion coefficient, and it is likely that 
Stokes-Einstein type models agree with this result for at least some 
classes of liquids. Henceforth, we will be concerned exclusively with 
the Free Volume Model, since the volume change mechanism allows us to 
compare magnetic and low gravity effects most directly. 

Volume Changes Due to Magnetic Fields and Low Gravity 

Let us begin by considering magnetic field effects. As usual, there are 
two ways to approach this problem, the first being through the theory 
of magnetostriction (Reference 21). This theory identifies three 
different components of magnetostricti ve volume change occuring over 
different ranges of field strength (see Figure 5). The simplest of 
these components, and the one occuring at lowest values of H, is called 
the form effect. In the form effect the fractional volume change is 
given by: 


M = 
V 


2 2 

NX^H^ 
JL 


2cy 


0 


( 88 ) 



FIGURE 5: MAGNETOSTRICTION EFFECTS 
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where N is a demagnetizing factor 
Xm is magnetic susceptibility 
G is the bulk modulus of elasticity, and 
H is the field strength. 

Although the other components are more complicated, it is obvious from the 
figure that ^V/V is always positive and generally increases with increas- 
ing H. 

The next approach to the problem of magnetic field-induced volume change 
is through the definition of isothermal compressibility, equation 12. 
Rewriting this equation in terms of finite volume and pressure differences, 
we have: 


i-^)j = 


(89) 


where 6 is the isothermal compressibility coefficient. It will be re- 
called that magnetic pressure was discussed in Section 2.2. The change 
in magnetic pressure for the situation of interest here may be expressed 
as (Reference 22): 

[Ho^ - (90) 

where Hq is the initial field strength before the liquid was placed in 
the field, and H is the resulting steady-state of field after the liquid 
is placed in the field. Thus: 

(^)t = (^) 3,0 - XmH"] (91) 

which is always positive for H < Hq. This will occur in cases where eddy 
currents or viscous flow are balanced so that the electric currents in 
the liquid are negligible. 

Now we consider the volume change resulting from removing a liquid quan- 
tity from the earth's surface and placing it in a region of low gravity. 
Once again we will use the compressibility equation, but first the 
pressure change must be found. Figure 6 depicts a container of liquid 
with a crystal being pulled from it. The crystal -melt system is 
initially on the earth's surface. 
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FIGURE 6: CRYSTAL-MELT ARRANGEMENT 


If the equilibrium pressure at the surface is Pg, the pressure at a point 
a distance z below the surface is given by (Reference 23): 

Pz = Ps + pgz • (92) 

If the liquid is now taken into space where there is no gravitational con- 
tribution, the difference in pressure at the example point will be: 

aP = ?l - P^ = P3 - (Ps + pgz) = -pgz . (93) 


Therefore, from equation 89: 

{^)j = 3pgz . (94) 

Since all terms in equation (94) are positive, it is evident that a 
decrease in hydrostatic pressure produces an increase in the volume of 
the liquid. This general result is of course confirmed by everyday 
experience with fluids at pressures above one atmosphere, and should be 
expected to apply between zero and one atmosphere as well. One must be 
careful, however, to interpret the factor g in equation (94) as the 
difference between the acceleration due to gravity at the earth's surface, 
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and the acceleration level in the space laboratory. This is due to the 
difficulty in obtaining absolutely acceleration-free environments even 
in an orbiting laboratory. But since acceleration levels are commonly 
on the order of 10"®g to 10" ^g in orbiting spacecraft, taking g to be the 
same as the earth surface gravitational acceleration is a good approxi- 
mation. 

Thus the application of a magnetic field or the removal of a gravitational 
(acceleration) field will cause a positive increase in liquid volume, 
aV. From the definition of free volume (equation 13) it is obvious that 
this volume increase is the same as an increase in the total free volume, 
A(Nvf). Therefore, according to. equation (87), both magnetic fields and 
low gravity conditions will cause an increase in diffusion coefficient, 
i.e., D' >0. 

Solidification Rate Dependence on Diffusion and Free Energy 

In the Free Volume Model, the crystal growth velocity or solidification 
rate of a crystal growing from an undercooled melt (Reference 2) is given 
by equation 28 which is repeated here: 


Uc 


!o 

R 


AG/kT 

e ] 


(28) 


For fixed fg and Rq, the now established fact that D increases in a 
magnetic field or low gravity conditions would indicate a corresponding 
increase in Ug if it can be shown that the exponential term decreases. 


For the magnetic case it has already been shown that T' > T. The free 
energy change is (Reference 24) : 


aG' = aG - It • i = a6 - Y 


(95) 


where c' is a constant. Since c', T and are positive, aG' < aG so that 


g AG'/kr < g AG/kT 


Linder low gravity conditions, there is no gravity field-dependent change 
in temperature to consider, only the change in free energy. This is: 


AG' = A6 + AVaP = AG - eV(AP)2 


(96) 
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where V is the volume at the earth's surface. Since (aP)^ is positive as 
well as B and V, obviously: 

AG' < AG ^ < g^G/kl 

for the gravitational field case. 


Hence the solidification rate, Uc» is increased by both the application of 
a magnetic field and by processing in a zero gravity field region. This 
is a purely formal and qualitative result, and a true comparison of 
solidification in space and solidification in magnetic fields must wait 
until the parameters required by equations 28, 87, 88, 91 and 95 can be 
obtained. However, it should be noted that a magnetic field-induced 
increase in crystal growth rate has already been observed in the 
laboratory (Reference 25), It is recommended that further research 
place emphasis on the quantitative comparison of solidification in space 
and solidification in the presence of magnetic fields, and on the com- 
parison of the increased growth rate predicted by equation 28 and that 
observed by Schieber (Reference 25). 
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APPENDIX A 

DERIVATION OF THE MAGNETIC FIELD EQUATION (49) 

The last two of Maxwell's equations of electrodynamics state that: 

V X 1 (A-1) 

? X f = -p ll (A-2) 

and Ohm's law for the case of a moving conducting liquid is: 

aP = J (A-3) 

where a is the electrical conductivity, the field f is referenced to the 
"fixed" frame of an observer and f' is the electric field moving with the 
l^iquid. Thus for a liquid of permeability y and moving with velocity 
u, we have: 


P = f + y u X ff . 


(A-4) 


Combining equations A-2, A-3 and A-4 yields: 

i V X J - y V X (u X It) 


^ 9H 

3t 


(A-5) 


or,^ubsti;^uting from A-1 and applying the vector identity v x v x H = 
-v^H when v • ft = 0 (as is the case in electromagnetism) , we find that: 


= -1 v 2 ]t + I X (u X ]t). 

9t ya 

If we now define the magnetic diffusivity to be: 
1 


V = 


ya 


and the magnetic Reynolds number to be: 
Rjyi — LUyo 


(A-6) 


(A-7) 


(A-8) 


where L is a characteristic length for the liquid system and u is the 
magnitude of u, equation A-6 becomes: 

which is equation 49. 


A-1 



